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Abstract 

The goal of this section is to assess the possibility that quarkonium production 
rates may be enhanced in nucleus-nucleus interactions at the LHC relative to 
that predicted by extrapolation of processes thought to be dominant at lower 
q ■ energy. This enhancement could follow from the effects of incoherent recom- 

CN \ bination mechanisms involving uncorrected pairs of heavy quarks and anti- 

' quarks which result from multiple pair production. Two different approaches 

^ ■ have been considered: statistical hadronization and kinetic formation. Updated 

predictions relevant to Pb+Pb collisions at the LHC are given. 

1 1. INTRODUCTION 

> : 

The utility of heavy quarkonium production rates in nuclear collisions as a signature of color deconfine- 
ment was proposed more than 15 years ago Q. Since one expects that the long-range color confining 

■ potential will be screened in a deconfmed medium, the quark and antiquark constituents of bound states 
will be liberated. As the system expands and cools, these constituents will, in general, diffuse away 
from each other to separations larger than typical hadronic dimensions. When the confining potential 
reappears, a given heavy quark will not be able to "find" its heavy antiquark partner and form heavy 
quarkonium. It must then bind with one of the antiquarks within range at hadronization. Since these an- 
tiquarks are predominantly the lighter u, d, and s flavors, the final hadronic states will preferentially be 
those with "open" heavy flavor. The result will be a decreased population of heavy quarkonium relative 

■ to that which would have formed if a region of deconfinement had not been present. This scenario as 
J> , applied to the charm sector is known as J / ip suppression. 

•i-H . 

^ . At LHC energy, perturbative QCD estimates predict that hundreds of pairs of charm-anticharm 

quarks will be produced in a central lead-lead collision. This situation provides a "loophole" in the 
Matsui-Satz argument since there will be copious numbers of heavy antiquarks in the interaction region 
with which any given heavy quark may combine. In order for this to happen, however, one must invoke 
a physical situation in which quarkonium states can be formed from all combinations of heavy quarks 
and antiquarks. This of course would be expected to be valid in the case that a space-time region of color 
deconfinement is present but it is not necessarily limited to this possibility. 

One can make a model-independent estimate of how such a "recombination" mechanism would 
depend on nuclear collision observables. For a given charm quark, the probability V to form a J/ip is 
proportional to the number of available anticharm quarks relative to the number of light antiquarks, 

JVg_ N c - C 

u,a,s ul 

In the second step, we have replaced the number of available anticharm quarks by the total number of 
pairs initially produced, assuming that the total number of bound states formed remains a small fraction 
of the total cc production. We normalize the number of light antiquarks by the number of produced 
charged hadrons. Since this probability is generally very small, one can simply multiply by the total 
number of charm quarks, N c , to obtain the number of J j ip expected in a given event, 

N - 2 

N m oc ^- , (2) 



where the use of the initial values N C c = N c = Nc is again justified by the relatively small number of 
bound states formed. 

The essential property of this result is that the growth of Nji^, quadratic in total charm, with en- 
ergy [ 2 ] is expected to be much faster than the growth of total particle production in heavy ion collisions 
0. Without this quadratic mechanism, J ftp production is typically some small energy-independent frac- 
tion of total initial charm production \4\. We thus anticipate that the quadratic formation will become 
dominant at sufficiently high energy. Generic estimates of the significance of this type of formation 
process can be made [5]. Here we look at specific predictions of two models which share the above 
properties and update the expectations to LHC energies. 



2. Statistical Hadronization 

The statistical hadronization model is motivated by the successful fits of relative abundances of light 
hadrons produced in high energy heavy ion interactions according to a hadron gas in chemical and 
thermal equilibrium [6]. Extension of the model to hadrons containing heavy quarks undeipredicts the 
observed abundances. This effect may be attributed to the long time scales associated with thermal pro- 
duction and annihilation of heavy quarks. The statistical hadronization model as first formulated for 
charm quarks |7| assumes that the cc pairs produced in the initial hadronic interactions survive until 
their subsequent hadronization, at which time they are distributed into hadrons according to the same 
thermal equilibrium parameters that fit the light hadron abundances. Chemical equilibrium abundances 
are adjusted by a factor 7 C which accounts for the non-thermal heavy quark density. One power of this 
factor multiplies a given thermal hadron population for each heavy quark or antiquark contained in the 
hadron. Thus the relative abundance of the J /if) to that of D mesons, for example, may be enhanced in 
this model. 

The value of 7 C is determined by conservation of the heavy quark flavor. For the charm sector, the 
conservation constraint relates the number of initially-produced c-cbar pairs N C £ to their distribution into 
open and hidden charm hadrons, 

N C c = "Terpen + Tc^hidden, (3) 

where N open is the number of hadrons containing one core quark and iVhidden is the number of hadrons 
containing a cc pair. For most applications, A^idden ( an d also multi-charm hadrons) can be neglected 
compared with N open due to the mass differences. Thus the charm enhancement factor is simply 



2N, 



1c = f^, (4) 

1 v opcn 

leading directly to the quadratic dependence of the hidden charm hadron population on N cc . One can 
then express the total number of J/ip in terms of the various thermal densities, n^, and the total number 
of cc pairs, One factor of system volume V remains implicit here. It is generally replaced by the 
ratio of number to density for total charged hadrons, n^/N^. Then the number of J ftp produced obeys 
the generic form anticipated in Eq. $2%. 

Nj = f* n J/* N cc (5) 
^open^ ^ch 

For collider experiments such as those at the LHC and RHIC, relating the corresponding central 
rapidity densities will be more relevant. Since Eq. © is homogeneous in the total particle and quark 
pair numbers, it will also be valid if these are replaced by their rapidity densities. To get an order of 
magnitude estimate, we choose a "standard" set of thermal parameters, T = 170 MeV and fiB ~ 0, for 
which the thermal density ratio is approximately 0.5. For a specific normalization, we assume dN^/dy 



= 2000 for a central collision at the LHC and take the initial charm rapidity density to be dN c z/dy = 25, 
roughly corresponding to N C c = 200 for central collisions (b = 0). Using these inputs, one predicts 
dNjj^/dy = 0.625, indicating that several J/ip will form through statistical hadronization in a central 
collision. To put this number in perspective, it is revealing to form the J/ip to rapidity density ratio, 
0.025 with the same assumptions. For comparison, one expects the corresponding hadronic production 
ratio to be of 0.01. This number would then be significantly reduced if placed in a region of color 
deconfinement. Thus the efficiency of J / ip formation via statistical hadronization at the LHC is expected 
to be substantial. 

These numbers can be easily adjusted to other charm and charged particle densities using Eq. 
Variations of the thermal parameters can also be investigated. For example, if the hadronization tempera- 
ture is decreased to 150 MeV, the prefactor combination of thermal densities increases by approximately 
a factor of two. 

The centrality dependence is controlled by the behavior of N C c and N^. The former should be 
proportional to the nuclear overlap function Taa{^) but is generally recast in terms of the dependence 
on the number of nucleon participants, -/V part . The calculation of -/V part requires a model calculation 
dependent on the total inelastic cross section, a in , as well as the nuclear geometry. We parameterize 
the expected behavior as a power-law N cc - oc iV part 4//3 . However, there will be deviations from this 
behavior for the larger values of a- m expected at the LHC [8]. The centrality dependence of at RHIC 
is also consistent with a power-law with exponent « 1.2 (3J[9). We will use the same dependence for our 
estimates at the LHC. It is clear that for sufficiently peripheral collisions one will encounter situations in 
which the average number of initially produced cc pairs is of order unity or less. At this point, one must 
revisit the assumptions of the original statistical hadronization model which assumed a grand canonical 
ensemble. The grand canonical approach is valid only when N C c is large enough for the fluctuations 
about the average value to be neglibile. Thus for peripheral collisions, one must recalculate the statistical 
results in the canonical approach where the charm number is exactly conserved, as noted in fTTH. Charm 
conservation can be implemented via a correction factor Iffl . 

N c ~ c = \ lc N op J^ cN J™\ + 7c 2 iV hidd en. (6) 

In the limit of large 7 c iV open , the ratio of Bessel functions approaches unity and the grand canonical 
result is recovered. In the opposite limit when 7 c iV open — > 0, the ratio of Bessel functions approaches 
57c-/Vo P en- In this limit, the dependence on N C c in Eq. (J5Jl changes from quadratic to linear. At the LHC 
this effect will not be relevant until one reaches very peripheral events, but at lower energies it can be 
significant over a much larger range of centralities fT2"ll . 

The results for dNj/^/dy as a function of iV part at the LHC are shown in Fig.^ The results are 
shown for three different values of dN C c/dy(b = 0), corresponding to iV C c(0) « 200, 150, and 100. 
There is a rapid increase with centrality due to the quadratic dependence of Nj/^ on N ce- 
lt is also interesting to look at these results normalized by iV par t, shown in Fig. [2] This ratio 
also increases with centrality, providing a signature for the statistical hadronization process that is less 
dependent on N cc - for the overall normalization. The corresponding results when normalized by dN C c/dy 
are shown in Fig. [3] The same general behavior is seen but the increase with centrality is less pronounced 
since the dN C c/dy is assumed to vary with a larger power, -/Vp art 4//3 . All of these ratios are at the percent 
level for central collisions and hence are larger than expected if the total J/ip population were due to 
initial production followed by suppression in a deconfined medium. 

The region of very peripheral collisions deserves some separate comments. First, there is a rise at 
low -/V part in both Figs. and [3] due to the onset of corrections from the canonical ensemble treatment. 
However, the extremely large values of the ratios as iVpart —> are an artifact of the decreasing inter- 
action volume, V — > 0. This calculation must be cut off before iV part = 2, i.e. only one interacting 
pair. It is also in this region where one must take into account the remaining J/ip from initial production. 
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Fig. 1: Statistical hadronization results for J/ip production as a function of Af par t at the LHC. 



Since the survival probability is maximum for very peripheral collisions and the statistical hadronization 
process is least effective in this same region, there will be a crossover in the relative importance of the 
two mechanisms. Some studies have already been performed for this situation at SPS and RHIC energies 

El HUE). 

Finally, there is another lower cutoff in centrality for the statistical hadronization results, needed 
to avoid a contradiction with the tpf/ (J/ip) ratio at the SPS. Since both of these states receive identical 
factors of y c , their ratio must be that predicted for chemical equilibrium in the absence of any charm 
enhancement or suppression. Although the measured ratio appears to be consistent for more central 
collisions [ 16 1, there is an indication that it rises sharply for more peripheral collisions. Most treatments 
have thus inserted a cutoff of iV part = 100, below which model predictions become inconsistent Q- 

The numerical values for dNj/^/dy are tabulated as a function of impact parameter in Table 1 for 
the three choices of initial charm multiplicity density and the default values of all other quantities. 



3. Kinetic Formation in a Deconfied Region 

The kinetic model has been developed [ 17 18 1 to investigate the possibility that J/ijj may form directly 
in a deconfined medium. This formation takes advantage of the mobility of the initially-produced charm 
quarks in a deconfined region. In order to motivate this view, consider the "standard" physical picture 
of deconfmement in which quarkonium is suppressed by collisions with free gluons in the medium [ 19 1. 
Then the formation process, in which a c and c in a relative color octet state are captured into a color- 
singlet quarkonium bound state and emit a color octet gluon, is simply the inverse of the the breakup 
reaction responsible for the suppression. This is an inevitable consequence of the suppression picture. 

The proper time evolution of the J/ip population is given by the rate equation 

dN J/rj} N c N-c 

—fa- = Ap 17(77 ~ XdNj /'^ Pa ' (7) 




Fig. 2: The statistical hadronization results divided by AT par t as a function of -/V part . 



where p g is the gluon number density and V(r) is the time-dependent volume of the deconfined spatial 
region. The reactivities Af,d are the reaction rates, (<7t> re i), averaged over the momentum distributions 
of the initial participants, i.e. c and c for Af and J/ip and g for Ad- 

The solution of Eq. (Q grows quadratically with N C c, as long as Nji^ <C N C c. In this case, we 

have 



Nj/^Tf) = e(r f ) Nj/^ro) +N%i' dr\ F [V(r) e(r)]- 

Jt 



(8) 



The function e(r^) = exp(— g?tAd p g ) would be the suppression factor if formation were neglected. 

The quadratic factor N^- is present, as expected, for the additional formation process. The normal- 
ization factor of iVch is not immediately evident, but is implicit in the system volume factor. This volume 
is now time-dependent, accounting for the decreasing charm quark density during expansion. Here the 
transverse area of the deconfined region is determined not just by the nuclear geometry but by the dy- 
namics which determine the extent of the deconfined region. This area is modeled by the energy density 
in terms of the local participant density in the transverse plane, n part (6, s = 0). The transverse area is 
denned by the ratio of the participant number to the local participant density. Note that the maximum 
local density is at s = 0. Thus, 

A T (b) = A T (0)[N pait (b)n paxt (0, s = 0)/iV part (0)n part (6, a = 0)] (9) 

These area effects will be more explicit when the centrality dependence is considered. 

The numerical results depend on a number of parameters, including the initial volume and temper- 
ature, the time expansion profile, the reaction cross sections, the behavior of the quarkonium masses and 
binding energies in the deconfined region, and the charm quark momentum distributions. For specifics, 
see Ref. [20|. Our previous results have used initial values N C c = 200, 150, and 100, spanning a rea- 
sonable range of expectations fTT\. The results are very sensitive to the initial charm quark momentum 




Fig. 3: The statistical hadronization results for J/ip production at the LHC, divided by the open charm multiplicity, dN^/dy, 
as a function of N pa , T t . 



distributions, as may be expected. We assume the charm p? distributions are gaussian and the charm 
rapidity distributions are flat over a plateau of variable width, Ay. The range 1 < Ay < 7 spans the 
range between an approximate thermal momentum distribution, Ay w 1, to a distribution similar to that 
of the initial pQCD production, Ay 7. 

The results as a function of the initial number of cc pairs produced in central collisions are shown 
in Fig. |4] There is a rapid decrease in formation with increasing Ay. The quadratic dependence on N C c 
is evident, but there is also a substantial linear component in some of the curves. This linear contribution 
arises because the final J/ tp formation by this mechanism is large enough for exact charm conservation 
to reduce the number of c and c quarks available to participate in the formation process. The curve 
labeled "Quadratic Extrapolation" uses a quadratic dependence derived from a fit valid only for low N C c. 
Note that the result for a thermal distribution is very similar to the assumption Ay = 1. 

The corresponding centrality dependence is presented in Fig. [5] where we give Nj^ at hadroniza- 
tion for three different initial charm quark momentum distributions, thermal, Ay = 4 and Ay = 7, as 
well as for our three choices of N C c(b = 0). 

Finally, the ratio of final J j ip to initial charm production is shown in Fig. |6] using the same 
parameters in Fig. |4] These ratios are most easily compared to either initial production or suppression. 
There is a substantial variation in the predictions and it is evident that a simultaneous measurement of 
open charm will be required for an interpretation. However, the centrality dependence is opposite to that 
expected in any pure suppression scenario. 

We have updated the calculations to include the charm quark momentum distribution from a lead- 
ing order pQCD calculation [22]. The rapidity distribution has a somewhat larger effective Ay and the 
Pt distribution does not fall as fast as a simple gaussian. As a result, the formation efficiency is further 
reduced. Such distributions may be most relevant, given preliminary results from RHIC l23ll24l . 
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Fig. 4: The J/ip production per central LHC collision in the kinetic model as a function of the initial number of cc pairs. 



The numerical values for Nj are compared with the statistical hadronization model results for 
dNj/^/dy in Tabled The overall magnitudes are comparable, although the centrality dependences 
differ somewhat. Thus details such as the resulting J/ip momentum distributions will be required to 
differentiate between these two models \22\. For completeness, Nji^ for the thermal distributions and 
the assumption Ay = 4 are presented in Table |2] 

4. Conclusions 

The "smoking gun" signature of the quarkonium formation mechanism is the quadratic dependence on 
total charm. For central collisions at the LHC one expects that this feature will lead to a total J/ij) rate 
greater than that produced by an incoherent superposition of the initial nucleon-nucleon collisions, even 
without any subsequent suppression due to deconfinement effects. In addition, the centrality dependence 
can be used to identify the quadratic dependence on charm assuming that the initial charm production 
scales with the number of binary collisions. Binary scaling leads to an increase of the ratio of J/ip to 
initial charm as the collision centrality increases, independent of specific parameters which control the 
overall magnitudes. A simultaneous measurement of total charm will be essential for such conclusions 
to be drawn. 

Uncertainties in the absolute magnitude of the formation process are inherent in the model param- 
eters. For statistical hadronization, one can constrain the thermal parameters to within a factor of two 
using the observed hadron populations. There is some additional uncertainty related to the lower cutoff 
on centrality needed to ensure the quarkonium ratios are consistent with an overall thermal picture. There 
is also the possibility that the correction for canonical ensemble effects will involve a thermal volume 
parameter not necessarily equal to the total system volume [25 j. In addition, the formation mechanism 
could be limited to those charm quarks whose phase space separation is within some maximum value, 
introducing another as yet unconstrained parameter [14|. With kinetic formation, a similar set of un- 
certainties exist. There are uncertainties in the space-time properties of the deconfinement region. In 
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Fig. 5: The centrality dependence of J/tp production in the kinetic model. 



addition, possible variations of charmonium binding energies and reaction cross sections in a deconfmed 
region are at present not well understood. There are indications that the efficiency of the formation 
mechanism is considerably reduced when included in a partonic transport calculation [26 1. 

The primary uncertainty in both models is still the initial number of charm quarks and their mo- 
mentum distributions. The tabulated J/ip results should be regarded in this light. Thus numbers may 
be only an order of magnitude estimate. However, the variation with centrality and total initial charm 
should provide experimental signatures which are largely independent of the overall magnitudes. 
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Table 1: Comparison of J/ip production at the LHC by the statistical hadronization (left-hand side) and kinetic formation 
(right-hand side) models. 
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Table 2: Kinetic J/ip formation at the LHC assuming both thermal charm momentum (left-hand side) and Ay — 4 (right-hand 
side). 





Nj/^ (Thermal) 
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